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Abstract 



We propose an interpretation of quantum separability based on a physical 
principle: local time reversal. It immediately leads to a simple characteriza- 
tion of separable quantum states that reproduces results known to hold for 
binary composite systems and which thereby is complete for low dimensions. 
We then describe a constructive algorithm for finding the canonical decom- 
position of separable and non separable mixed states of dimensions 2x2 and 
2x3. 
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Typeset using REVTgX 



Entanglement, inseparability and nonlocality are some of the most genuine quantum 
concepts. As it has been pointed out, we still lack a complete classification of quantum states 
in local and non-local ones , or more precisely we still lack the complete understanding of 
non-locality. Thus, while for pure states it is well established since long ago that the non-local 
character of the composite system is revealed in different but equivalent ways, the situation is 
drastically different for mixed states. For example, for pure states the violation of some kind 
of Bell inequalities [0 , or the demonstration that no local hidden variable models can account 
for the correlations between the observables in each subsystem, are equivalent definitions 
of non-locality. But for mixed states, described by density matrices, such equivalences fade 
away. Consider a composite quantum system described by a density matrix p in the Hilbert 
space Tia ® 'Hh- In the frame set by the concepts of our starting sentence, product or 
factorizable states are the simplest possible. They are of the form Pp = Pa® Pb, i-e. for 
them, and only for them, the description of the two isolated subsystems is equivalent to 
the description of the composite system. Recalling that subsystems are described by the 
reduced density matrices obtained via partial tracing: pa = Tri,p ( p^ = Tvap), a density 
matrix corresponds to a product or factorizable state if and only if 

p = Tr-bp Trap P = Pp (1) 

Also their index of correlation defined in terms of von Neumann entropies of the system and 
subsystems, 

Ic = Tr p In p-Tr Pain Pa -Trpblnpb (2) 

vanishes, and this happens only for them . Their subsystems are uncorrelated. Any state 
which is not a product state presents some kind of correlation. They are called correlated 
states. Quantum mechanics has taught us that there is a hierarchy of correlations, and the 
physics in the different ranks is different. The simplest correlated systems are the classically 
correlated or separable systems. Their density matrices can always be written in the form: 

Ps = ^PiPai ® pu] 1 > Pi > 0; = 1 (3) 

i i 

i.e. as a mixture of product states. Their characterization is notoriously difficult. Thus, 
given a density matrix which is known to describe a separable system no algorithm for 
decomposing it according to eq. is known; besides, the decomposition is not unique. In 
fact, only very recently Peres and the Horodecki family |^,^ have obtained a mathematical 
characterization of these states, at least when the dimension of the composite Hilbert space 
is 2 X 2 or 2 X 3. For these cases the necessary and sufficient condition for separability is that 
the matrix obtained by partially transposing the density matrix p is still a density matrix, 
i.e. hermitian, with unity trace and non-negative eigenvalues 

P^^ = (P^")*>0 ^ p = ps (4) 

For composite systems described by Hilbert spaces of higher dimensions, the positivity 
condition of p^^ is only a necessary one for separability ||^. Following the hierarchy of 
correlations, we find the states that are no longer separable (or classically correlated ), 
i.e. p 7^ Ps- These states are called "EPR-states" 0, "inseparable", "non-local", and 
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sometimes "entangled" or simply "quantum-correlated" to emphasize that their correlations 
are not strictly classical anymore, though often these labels do not refer to exactly the same 
states. This confusion reflects the need of a further subclassification of the inseparable states 
according to whether they admit local hidden variables, whether they violate some kind of 
Bell inequality etc.. 

The aim of this Letter is threefold. Firstly, we give a physical interpretation of the 
mathematical characterization of separability (Eq. (^). Secondly, we provide a constructive 
"canonical" algorithm for decomposing any separable matrix (of dimension < 6 ) into a very 
small finite incoherent sum of product vectors. Finally, we show that a similar decomposition 
holds for inseparable states, with the signature of non-separability being expressed by some 
non-positive weights in the canonical decomposition. 

Let us first analyze the problem of separability from a physical point of view. We start 
by considering symmetry transformations in the Hilbert space of each subsystem. We limit 
ourselves in this paper to just binary composite systems , i.e. Ti = Ha ® Hb- Wigner's 
theorem tells us that every symmetry transformation should always be implemented by a 
unitary [U) or antiunitary {A) matrix. The direct product of unitary matrices Ua Ub 
(or antiunitary matrices Aa Ab) is a unitary (or antiunitary) matrix in the Hilbert space 
of the composite system 7i, and one can give an unambiguous definition of how such a 
transformation acts on any ket G 7i. However, the combination of a unitary and 
an antiunitary transformation Ua (S> Ab (or Aa ® Ub) results in a transformation which is 
neither unitary nor antiunitary in Ti, whose action on a general ket of the composite system 
1^) e H, furthermore, cannot be properly defined. However, its action on a product ket 
|e) ® I/) = |e, /), (where |e) G Ha and |/) G Hb) is, but for a phase ambiguity, well defined. 
Thus, the action of combined transformation of the type Ua®Ab on projectors corresponding 
to pure product state is well defined without any ambiguity. As a separable state can always 
be rewritten as a statistical mixture of product vectors: 

i i 

it is clear that under the combined transformation Ua ® Ab (or Aa Ub) Ps transforms into: 

Ps-^Ps=j:p^{\e[){e:\®\f:){f:\) (6) 

i 

where \e[) = Ua\ei) G Ha ; \f'i) = ^blfi) £ Hb- Therefore, p'^ describes also a physical 
state so that p^ is a positive defined hermitian matrix (with normalized trace). This is what 
characterizes separable states: that any local symmetry transformation, which obviously 
transforms local physical states into local physical states, also transforms the global physical 
state into another physical state. (Here and in what follows "local" means that it refers to 
the subsystems). 

There exists only one independent antiunitary symmetry and its physical meaning is well 
known: time reversal. Any other antiunitary transformation can be expressed in terms of 
time reversal (as the product of a unitary matrix times time reversal). We are thus proposing 
that quantum separability of composite systems implies the lack of correlation between 
the time arrows of their subsystems. In other words: for separable states, the state one 
obtains by reversing time in one of its subsystems is also a physical state. Loosely speaking, 
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systems which are classically correlated (separable) do not have memory of a unique time 
direction in the sense "EPR" correlated states have, and they are thus compatible with a 
time evolution which factorizes into the product of two opposed time evolutions. Changing 
the time direction in only one of the subsystems but not in both leads to a physical state 
since their time arrows are uncorrelated. 

We can define now the "separability" operator as the simplest possible transformation 
of this type: 

S = Ia(S)Kb (7) 

where la stands for the identity acting in the first subsystem Tia and Kb is the complex 
conjugation operator acting in Tib 0. It is straightforward to check that in the Hilbert- 



Schmidt basis [|T0| of 2 x 2 composite systems: 

SpS = p^" (8) 

for any p, whether it is separable or not, in spite of the fact that the action of 5 on a general 
I^E') G 7i cannot be properly defined. (This feature of being able to define a transformation 
on density matrices which one cannot define on kets is known for some nonunitary trans- 
formations, as e.g. a decohering time-evolution). Finally, as local time reversal is locally 
unitarily equivalent to expression (|^), eqs.(H) and (||) state that for 2x2 composite systems 
a state is separable if and only if it does not correlate local time flows. The same holds for 
a 2x3 composite system. 

Let us go to the second point of our Letter, the "canonical" decomposition of a separable 
density matrix. Until very recently it was not known whether one could always find a finite 
incoherent sum of pure product states for any separable ps- P. Horodecki |ll| and Vedral 
et al. |]T2[ have shown that any separable state ps can be decomposed into an incoherent 



sum of at most = {dim{Ha) x di'm(Hb)Y pure product states, although no algorithm for 
obtaining this decomposition seems to be known. We will limit ourselves here, again, to the 
simplest possible case, i.e. binary composite systems of dimensions dim(Ha) = dim(Hb) = 2 
. In such cases, any separable density matrix can be written as a statistical mixture of at 
most N=1Q pure product states. Let us show here that one can do much better: indeed, 
any separable density matrix can be written as a convex combination of at most = 5 
pure product vectors. The whole proof of such a decomposition is based on the following 
theorems: 

Theoreml For any plane Vi in ® defined by two product states |f i) and \v2) (where 
= |ej) (8> \fi); |ej) G Ha and G Hb) either all the states in this plane are product 
states, or there is no other product states in it. 

Theorem2 There exist planes V2 in ® which contain only one product state. 
Theorem 3 Any plane P3 in ® contains at least one product state. 
The proofs of the theorems are simple. It is convenient to express, with the help of the 
SU{2) ® SU{2) transformations, the planes defined by the theorems (denoted as Vi,V2 and 
V3) as: 
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with < A,C <7i/2; < B,D <27i, and ai, [3i G C. 

I 



,0. 



.0. 





cosA 



\ 



e**" sin A cos B 



\ sin A sin 5 / 

with < A < 7r/2; < 5 < 7r/2; < C < 27r, and ^2, e C. Finally: 



^ cos v4 



V sin A y 



/ sin A cos -B \ 
sin B cos C 
e*^ sin 5 sin C 
V — cos A cos B ) 



(10) 



(11) 



with < A < 7r/2; < 5, C < 7r/2; < D < 27r, and as, e C. There are further condi- 
tions that have to be imposed in Eq.(|TT]) to ensure that the second state is not a product 
state. Solving for the values of a, and /3j allows to prove the theorems. A consequence of 
the above theorems is the following corollary: 



Corollary If p has rank 2 and is separable it can always be written as a statistical mixture 
of two pure product states, and thus p^^ is also of rank 2. 

Consider now a separable state pg such that both itself and its partially transposed matrix 
are of rank 4, 

APs) = r(pr^) = 4 (12) 
(All other cases are subcases of this one, as we shall see immediately). Define now 

p{p) = Y^{Ps-p\ei,h){eufi\); 0<p<l (13) 

where |ei) G Ha and G Hb are completely arbitrary kets. For p small enough both p 

and p"^* 

p'^''ip) = Yz^ip^'-p\^^^fi)(^^'fi\y' o<p<i (14) 

are positive, and therefore, due to eq.(||), separable. Let us denote by pi the smallest value 
for which a zero eigenvalue appears in p{p) or p^''{p). Let us assume that for pi one eigenvalue 
of p{p) is equal to zero, i.e. r{p{pi)) = 3 and r{p^^{j)i)) = 4 (the same argument holds for 
the opposite case). Consider now a new pure product vector state belonging to the range of 
p(pi); |e2, /2) £ '^(p(Pi)) and define a new density matrix: 

P(P) = T^(p(Pi) -P|e2,/2)(e2,/2|); 0<p<l. (15) 
I — p 

As before, for p small enough, both p{p) and p^*'{p) are non- negative and thus separable. 
Let us denote by p2 the smallest value of p for which either p{p) or p^''{p) develop a new 
vanishing eigenvalue. It cannot be p(p) unless, because of the corollary, p^*'{p) develops 
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simultaneously two vanishing eigenvalues. Therefore, it is in general p^*'{p) which will develop 
a new vanishing eigenvalue, so that 

^(p(P2))=r(/''(p2)) = 3. (16) 

As p(p2) has a decomposition of the type of eq.@ with at least three terms, and p^*'{j)2) 
has the corresponding partially transposed one, there exist always a product state satisfying 



TT| : |e3, /s) G lZ{p{p2)) and |e3, f^) G 'R-ip^'' {P2)) where the following identity: 

(i0)(0ir = (i0)(0ir (17) 

and some results by Hughston et al. [|13| have been used. Define now: 

P{P) = Y^{P{P2)-P\e3j3){e3,f3\y, 0<p<l. (18) 
It is clear from the above results that a p^ exists such that: 

r(p(p3) > 0) = r{p^''iP3) > 0) = 2 (19) 
Finally, from the corollary, we know that there exists always the following decomposition: 

P{P3) =P4|e4,/4)(e4,/4| 

+ (l-p4)|e5,/5)(e5,/5|; 0<p4<l. (20) 
And from there our final results follows immediately: 

Ps = PiPi+P2il -Pi)P2 
+ Ps{l-p2){l-pi)P3 

+ P4(l-P3)(l-P2)(l-Pl)P4 

+ (1 - p,){l - P3)(l - P2)(l - Pl)P5 (21) 

where Pi = |ej, fi){ei, fi\ are projectors onto pure product states. 

It is obvious that this decomposition is far from unique. It is also clear that often one 
can simplify the first two steps to one, i.e. go from expression ([l^) to expression ([TBD by 
subtracting only one single pure product state. This can always be done, by continuity, when 
for some |e, /) it happens that r{p{pi)) = 3 and r{p^''{pi)) = 4 and for others r{p{pi)) = 4 
and r{p^^{pi)) = 3. We do not know, right now, whether it can always be done, i.e. if any 
separable density matrix can always be written as a statistical mixture of just A^=4 pure 
product states. 

Let us now discuss our third result. Recall that for dim{l-C) = 4 a state has quantum 
correlations iff p^^ has at least one negative eigenvalue: 

infa{p^^)<{)^p = p, (22) 

where cr(p) means the spectrum of p and the subscript q means quantum correlated. Let us 
prove now that there is only one negative eigenvalue. If there were two, one could always find. 
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according to our theorem, a product state |e, /) in the plane defined by the corresponding 
two eigenvectors, and for which obviously 

{e,f\p^''\ej)<0. (23) 
But the above expression is equivalent, recalling the Hilbert-Schmidt decomposition, to 

(e,/1p,|e,n <0, (24) 
which is impossible. Let us then define 



1 

I+P1+P2 

+ P2|e2,/2)(e2,/2|),P^>0, (25) 



P^'iPuP2) = , , ^ , ^ (P^" +Pi|ei,/i)(ei,/: 



where |ej, /j) are the two product states of the Schmidt decomposition of the negative eigen- 
value eigenvector of p^*"; \n) = Ci|ei, /i) + C2|e2, /2)- For some finite values of pi and p2, pi 
and p2, r{pg''{pi,p2) > 0) = 3. As Pq{pi,P2) > the algorithm proceeds as before for the 
separable states. Thus 

Pg = (1 +Pi +P2)p.(4) -pi|ei,/*)(ei,/*| 

-P2|e2,/2^)(e2,/2l, (26) 

where Ps(4) is a statistical mixture of = 4 pure product states. Certainly, expression p6| ) is 
not a statistical mixture in the sense that two weights in the decomposition are negative, but 
that only means that Pg, which is a statistical mixture of pure states, is inseparable. Often 
one can find a decomposition of the type ( PB| ) but with only five terms, either having Ps(3) 
or only one p, but we do not know yet if this is always possible. Also, pi are measures of the 
inseparability of the state, which, supplemented with an adequate minimization procedure, 
might lead to a faithful quantification of entanglement p^ , p3| . 

Finally, let us illustrate our procedure with a simple example. Consider a pair of spin-| 
particles in an impure state consisting of a singlet fraction x and an isotropical mixture of 
the singlet and the triplet mixed in equal proportions 0: 

p^ = x|vI/-)(M/-| + il^(|M/-)(M/-| 

+ |^+)(^+| + |$+)($+| + |$")(<^"|) (27) 

where < x < 1, = 1/^2(1 ti)±| it)), and = l/v^(| TT)±| ii))- The condition 
of separability (eq. ^ shows that p^ is separable for a; < 1/3 and inseparable otherwise. 
One simple decomposition by following the procedure we have indicated is given by: 



_ {l+3x)(l-x) 
Pi — 4{l+x) 

(l-3x)(l+x)^ 
P2 - — 



(3+2x+3x2)(l-a;) 



P3 = i 

P4 = l 



|ei,/i) 



i^3,/3) = ^(_^)^^(:3i)^ (28) 

|e5,/5) = \e4.y <S) 1/4)* 
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which also holds when is inseparable, in which case p2 becomes negative. One can extend 
most of our results straightforwardly to dim\H]=Q, but not beyond, as Eq.(||) does not hold 
anymore. 

To summarize, we propose to characterize separable states by their inability of correlating 
local time arrows. For low enough dimensions this characterization is complete. We have 
also shown that it is always possible to find a decomposition with at most five pure product 
states to express any separable density matrix (for dim\H] = 4). Moreover, when the state 
is inseparable a similar decomposition with at most six pure product states holds, where 
one or two of them have now a negative weight. We believe our results are a step forward 
in the understanding of quantum (non)separability. 
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